Disorder- quenched Kondo effect in mesosocopic electronic systems 
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Nonmagnetic disorder is shown to quench the screening of magnetic moments in metals, the Kondo effect. 
The probability that a magnetic moment remains free down to zero temperature is found to increase with dis- 
order strength. Experimental consequences for disordered metals are studied. In particular, it is shown that 
the presence of magnetic impurities with a small Kondo temperature enhances the electron's dephasing rate at 
low temperatures in comparison to the clean metal case. It is furthermore proven that the width of the distri- 
bution of Kondo temperatures remains finite in the thermodynamic (infinite volume) limit due to wave function 
correlations within an energy interval of order 1 /r, where r is the elastic scattering time. When time-reversal 
symmetry is broken either by applying a magnetic field or by increasing the concentration of magnetic impuri- 
ties, the distribution of Kondo temperatures becomes narrower. 
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I. INTRODUCTION 

A local magnetic impurity changes the ground state of a 
Fermi liquid due to the correlations created by the antifer- 
romagnetic exchange interaction between its localized spin 
and the delocalized electrons. 1,2 At temperatures below the 
Kondo temperature Tk the spin of the magnetic impurity is 
screened by the formation of a singlet state with the conduc- 
tion electrons. 3,4 Disorder affects the formation of this Kondo 
singlet in various ways. 5 In particular, the Kondo temperature 
may depend on the positioning of the magnetic impurities in 
the host lattice. Thus, for a sample containing many mag- 
netic impurities, Tk may be distributed due to fluctuations of 
the microscopic exchange coupling. 6,7,8 Impurities (magnetic 
or nonmagnetic) and lattice defects also cause fluctuations in 
the local density of states (LDOS) of the conduction electrons 
at the magnetic impurity site. Thus, the distributions of the 
Kondo temperature and of the LDOS are related in a disor- 
dered metal .-2^ 

In this paper we explore the consequences of nonmagnetic 
disorder for the Kondo effect. We base our study on the sta- 
tistical properties of Tk for the different dynamical regimes 
of a disordered metal. Note that the Kondo temperature ac- 
counts for a crossover rather than a sharp transition in the 
metal properties and that may raise the question whether Tk 
is sufficiently well defined. In weakly disordered metals it is 
possible in principle to extract Tk by fitting physical quanti- 
ties such as the spin susceptibility and the electron dephasing 
rate versus temperature against the corresponding universal 
scaling function. While that scaling function may be modified 
itself by the disorder, 8,11 the Kondo screening is still governed 
by Tk, and one may study its sample-to-sample fluctuations. 
Below, we show that these fluctuations may have a significant 
impact on the low-temperature properties of a metallic sample 
even in the thermodynamic, infinite volume limit. 

It is useful to compare the magnitude of Tk to other rel- 
evant energy scales of a metallic sample. These scales are 
the mean level spacing A = l/vL d , the Thouless energy 
E c = D e /L 2 , and the elastic scattering rate 1/r. Here, v 



denotes the density of states at the Fermi level, D e 



?r/d 



is the diffusion constant, and d is the dimensionality of diffu- 
sion (we have set h = ks = 1). This comparison allows one 
to establish several distinct regimes, as sketched in Fig. [TJ We 
note that, in practice, Tk does not exceed 1 /r in metals, there- 
fore standard perturbation theory in the disorder potential may 
not be used to describe the effect of disorder on the Kondo 
effect. However, there is a large regime of experimental in- 
terest where a diagrammatic expansion may be used, namely, 
E c < T K < l/r^For small samples, when A < T K < E c , 
random matrix theory (RMT) can be applied and the distribu- 
tion of Kondo temperatures is expected to scale with A alone. 
In the opposite limit, when the sample is so large that the lo- 
calization length £ is smaller than the linear size L, the Thou- 
less energy and the global level spacing are irrelevant and the 
Kondo temperature is determined by the mean level spacing 
for states localized in the vicinity of the magnetic impurity, 
A^ = 1/V£ d i 13,14 ' 15 In this case, one can expect the distri- 
bution of Kondo temperatures to scale with A^ instead of A. 
Finally, when Tk is smaller than the spacing between neigh- 
boring energy levels at the Fermi surface (either A or A/), 
the distribution of Tk will be mainly determined by the fluc- 
tuations of the wave functions and level spacing of the two 
eigenstates closest to the Fermi energy. Accordingly, in this 
regime the Kondo temperature is determined by the coupling 
of the magnetic impurity to a two-level system (2LS). 

Our analytical calculations show that the average Kondo 
temperature is enhanced by weak disorder. Disorder also in- 
duces a finite width in the distribution of Tk which is shown 
to persist even in the infinite sample, thermodynamic limit. 
This effect is due to the existence of local correlations between 
eigenfunctions at different energies within a macroscopically 
large energy interval 1/r. These spectral correlations also 
lead to a bimodal structure in the distribution of Tk, with a 
long tail stretching toward Tk = 0. The numerical simula- 
tions show that for a sizeable number of disorder configura- 
tions, no finite value of Tk exist within the lowest order self- 
consistent perturbation theory, indicating that in those cases 
the magnetic moment is not screened by the conduction elec- 
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FIG. 1: Overview of the different regimes of a mesoscopic Kondo 
system of linear size L, as parametrized by the system's mean level 
spacing A = l/vL d and the unaveraged Kondo temperature of a 
single magnetic impurity Tk- E c is the Thouless energy, 1/r is the 
elastic scattering rate, and = , with £ being the localiza- 

tion length. The acronyms denote: Two levels coupled to the mag- 
netic impurity (2LS), diagrammatic expansion in diffusion ladder di- 
agrams (DD), and perturbative regime (PT), with 1/g — (Eft) -1 
as the small expansion parameter. We only consider metallic systems 
where the relation Aj < 1/t < Ef is fulfilled. 



irons. We find that there is a finite probability of having free 
moments even for values of J such that the average Kondo 
temperature (Tk) is much larger than the mean level spacing 
A. 

For small metallic grains, it is well-know that RMT pro- 
vides a very good description of the statistical fluctuations of 
single-particle states. We therefore use RMT to investigate 
analytically and numerically the statistics of Tk in such sys- 
tems. The ratio between the standard deviations of Tk for 
the orthogonal and unitary ensemble is found numerically to 
be equal to \pl for a wide range of exchange couplings J, 
in agreement with analytical calculations. Furthermore, it is 
shown numerically that a finite number of free moments per- 
sists to exist in this regime. This result is consistent with an 
analytical calculation based on the 2LS approximation. 

In regard to the width of the distribution of Tk in the 
RMT regime, we find a counter-intuitive result at first glance, 
namely, that the distribution becomes narrower when increas- 
ing the concentration. This effect maybe understood as a 
consequence of the suppression of the weak localization cor- 
rection by the breaking of time-reversal invariance, which in 
RMT is seen as a reduction in wave function fluctuations. We 
also consider whether additional correlations between wave 
functions induced by breaking of time-reversal symmetry 16 
within an interval of order of the spin scattering rate l/r s can 
lead to an enhancement of the width of the distribution of Tk- 
We find that while the width increases by a factor \/Tkt s as 
compared to the pure RMT ensembles, it nevertheless van- 
ishes in the thermodynamic limit when A — ► 0, and therefore 
is negligible in comparison with the magnitude of the broad- 
ening caused by wave functions correlations beyond RMT. 

The large fluctuations of Tk in disordered metals impact 
the behavior of several thermodynamic and transport proper- 
ties at low temperatures. We consider in particular the influ- 



ence of these fluctuations on the temperature dependence of 
the electron's dephasing rate. We find that nonmagnetic disor- 
der combined with spin scattering due to magnetic impurities 
increases the dephasing rate at temperatures T < (Tk}- The 
dephasing rate is readily accessible in transport experiments 
by measuring weak-localization corrections to the resistance. 
Our predictions are in qualitative agreement with recent ex- 
perimental results. 

Our work details and extends the discussions provided in 
a previous letter (Ref. 1 1 CM) and also includes new numerical 
and analytical results. The presentation is organized in the 
following way. In Sec. [II] we introduce the problem of find- 
ing the distribution of the Kondo temperature as defined by 
an integral equation. The central analytical and numerical re- 
sults of our paper are presented in Sees. UTTI and [TVl respec- 
tively. There, the distribution, average, and standard deviation 
of the Kondo temperature are studied as function of disorder 
strength by taking eigenstates and eigenenergies from a disor- 
dered tight-binding model. The average and standard devia- 
tion are derived analytically as function of A, E c , and the dis- 
order strength and compared with the numerical results. We 
also study how the number of free moments depends on the 
exchange coupling and disorder strength. 

In Sec. IVII we consider the RMT regime. We present the 
analytical results obtained for the distribution of Tk by taking 
into account only wave function fluctuations. We compare an- 
alytical estimates of the average and standard deviation of Tk 
with numerical results obtained from simulations of random 
matrix ensembles. Both large- and small- Tk limit are consid- 
ered. 

The dependence of the distribution of the Kondo temper- 
ature on the concentration of magnetic impurities is studied 
in Sec. I VIII In Sec. IVIIII we describe how the fluctuations 
of Tk affect the dephasing time of electrons in metallic thin 
films and quasi-one-dimensional wires. Finally, in Sec. |TX] 
we summarize our results and make our concluding remarks. 
Appendix A includes the derivation of the integral equation 
for the Kondo temperature in a one-loop approximation for a 
disordered sample, together with a critical discussion of its va- 
lidity. In Appendix B, for reference, we briefly present the so- 
lution of the integral equation in the clean limit using a picket- 
fence spectrum (equidistant energy levels). 

II. FORMULATION OF THE PROBLEM 

We will consider a single magnetic impurity in an isolated, 
phase-coherent metallic sample where the energy levels are 
well separated, a system previously referred in the literature 
as the Kondo boxji2J£ To determine Tk, we will adopt the 
following self-consistent equation j 10 ' 19 i 20 

l = lf fl '^'' tanhfg^V (1) 
2N ^E n -E F \ 2T K J 

where N is the number of states in the spectrum, ft = L d is 
the sample volume, Ep is the Fermi energy, and J is the ex- 
change coupling between the magnetic impurity and the delo- 
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calized electrons. The eigenenergies and eigenfunctions of the 
sample are, respectively, E n and ip n (r), while the impurity is 
located at r = 0. Equation (Q~|) is obtained from second-order 
perturbation theory in Appendix lAl A similar expression can 
be obtained from the zero-temperature, self-consistent solu- 
tion of the one-loop renormalization group (RG) equation.^ 2 - 1 - 
However, in that case one misses the tanh factor (which ac- 
counts for the finite-temperature occupation numbers) and the 
result is only valid for Tk 3> A. While the approximations 
involved in deriving Eq. ([TJ are not sufficient for describ- 
ing the properties of the system below Tk, it is important 
to remark that Eq. (Q]l does yield a good estimate for the 
Kondo temperature, which is the relevant scale for the low- 
temperature behavior of Kondo systems. The two-loop cor- 
rection has been found to change the Kondo temperature by 
a factor y/ J/D& In the thermodynamic limit (infinite vol- 
ume), the physics at temperatures T <C Tk is known to be 
that of an effective Fermi liquid, where the Kondo temper- 
ature determines the Landau parameters. 3 For example, the 
effective mass, the density of states, and thereby the spe- 
cific heat become enhanced by a factor 1 + n m / (itvTk), 
while the paramagnetic susceptibility is enhanced by a factor 
1 + 2n m / {itvTk ), where n m is the concentration of magnetic 
impurities 

It is convenient to rescale Eq. ((TJ to a dimensionless form, 
in which case it becomes 



1 = 



1 N 

2x ^ s n \2kJ 



(2) 



where x = D/J, n = Tk/D, x n = f2|^> n (r)| is the proba- 
bility density of the nth eigenstate, and s n = (E n — Ep)/A is 
the corresponding eigenenergy measured relative to the Fermi 
energy in units of A. In the following, we will assume that 
the Fermi energy is in the middle of two energy levels, so that 
the number of electrons in the Kondo box is even. (For an odd 
number of electrons, the unpaired electron at the Fermi en- 
ergy forms a singlet with the magnetic impurity with binding 
energy J and the problem becomes trivial.) Using Eq. (O, the 
distribution of Kondo temperatures is determined by 



P{Tk) = 



dT K 



D/J 



dx (5(x-F(T K ))) 



(3) 



where 



1 N 

F(T K ) = -J2-^ 
I *■ — ' s 

n=l 



■ ' r. 



(.2k) 



We note that a similar expression appears in the calculation of 
the distribution of NMR- Knight shifts for metallic grains with 
a finite level width ri 22 i 23 ' 24 There, an exact solution of the 
problem was derived within RMT in the limit N 3> 1. How- 
ever, two aspects make the problem defined by Eq. (01 much 
harder to solve. First, in the NMR-Knight shift case, instead 
of Eq. the calculation involves a sum over terms which 
behave as l/(s 2 + T 2 ), decaying faster than l/s n at large 



energies. Second and more importantly, since the Kondo tem- 
perature enters in Eq. © essentially as the low-energy cut- 
off of the sum, we cannot apply the techniques used in Refs. 
l22ll23U24l to the calculation of the distribution of Kondo tem- 
peratures. Thus, the derivation of an exact, closed expression 
for the P(Tk) is very nontrivial, even in RMT. 

In RMT, for (Tk) A, one would expect the distribution 
of Kondo temperatures to be close to a Gaussian. The reason- 
ing behind that is the following. In this limit one can go back 
to Eq. ©, expand the tanh factor and find that Tk is given 
approximately by a sum of random variables (the wave func- 
tion amplitudes). The central-limit theorem then ensures that 
a Gaussian distributed variable results from summing over in- 
dependent random variables. However, as we show below, 
even for large (Tk) one finds a non-Gaussian behavior at the 
tails of the distribution. This feature is enhanced when con- 
sidering disorder effects beyond RMT due to the appearance 
of local correlations between wave functions at different en- 
ergies. It is therefore clear that in order to obtain the correct 
distribution of the Kondo temperature in disordered metals, it 
is crucial to abandon the assumption of independent random 
wave functions in Eq. (f2]). The main effect of wave function 
correlations is to broaden the distribution of Tk- This broad- 
ening survives the large-volume limit. In Sec. |IV] we show 
that correlations also induce a bimodal structure in the distri- 
bution. This should be compared with the rather featureless 
distribution obtained when one uses uncorrelated RMT wave 
functions amplitudes to evaluate Tk 3 

For reference, the behavior of the Kondo temperature in the 
clean limit is discussed in Appendix iBl 



III. FLUCTUATIONS OF T K IN A DISORDERED METAL 

The effect of nonmagnetic disorder on the statistical fluctu- 
ations of the Kondo temperature is a difficult problem to study 
analytically. However, in the weak-disorder limit it is possible 
to derive expressions for the average and standard deviations 
of Tk beyond RMT by taking explicitly into account wave 
function correlations. To this end, we recall that the LDOS is 
defined as 



N 

E 

71=1 



p(r,£) = V|i(r)| 2 i(£-£„) 



(4) 



In our calculations, we will make use of the well-known cor- 
relation function of LDOS, as defined by 

R 2 (r,Lu) = (p(r,E)p(T,E + u)) 

= ±J2(\i>n(*)\ 2 \i>m(r)\ 2 d(E-E n ) 

n.m 



xS(E + u- E m )) 



(5) 



It is convenient to separate diagonal from off-diagonal terms 
in Eq. (|5). The correlation function of LDOS can then be 
rewritten in terms of the spectral correlation function ^(w), 

J2a(r,w) = R 2 (uj)n 2 (\Mr)\ 2 Hm(r)\ 2 ) U=E„-E m 

+ ^/A)fi 2 (|V„(r)| 4 ). (6) 
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Let us introduce the standard dimensionless parameter g = 
EpT to quantify the disorder strength. For u < E c , the func- 
tion i?2(<^) is an oscillatory decaying function with correc- 
tions of order 1 j g 2 to the leading RMT term. For the time- 
reversal symmetric case (J3 = 1), it reads 



1 



sin 2 (7rs) 



-sgn(s) - Si(7rs) 

0(i/g 2 ), 



cos(7rs) sin(7rs) 

7TS (ns) 2 



(7) 



where, Si(z) = L dy sin y/y.— For frequencies exceeding 
the Thouless energy, u> > E c , the oscillatory part of the spec- 
tral correlation function decays exponentially, while there is a 
correction of order l/g 2 which decays as 1 / s without oscilla- 
tions. While for pure RMT there are no correlations between 
wave functions at different energies, in systems with white- 
noise disorder these correlations are of order l/g, namely, 25 



n a (|V>„(r)| 2 |iMr)|> t 
for n ^ m, while 

n 2 (|^„(r)| 4 > = 



=E n -E„ 



1+-Ren(w) (8) 

P 



1 + -Re 11(0) 

P 



(9) 



These analytical results have recently been confirmed numer- 
ically in Ref. Note that the correlation is stronger when 
time-reversal symmetry is present (/3 = 1). The dependence 
on the disorder enters through the summation over diffuson 
modes, which is here represented by 



n(w) 



A 



7r ^— ' D e q 2 — ILO 
q 



In two dimensions (2D) and for L > I, one obtains^ 5 - 



n(^) M = J_ ln (i/i 2 -i»/D e 



2irg \l/L 2 -iuj/D t 



(10) 



(11) 



with I = vpT and D e — vfI/2. 

Since the Kondo temperature is defined by a sum over all 
eigenstates in the band, fluctuations of wave function ampli- 
tudes at the position of the magnetic impurity add up quite ef- 
fectively when the wave functions are correlated over a large 
energy range. These correlations are nonzero over an interval 
of the order of the elastic scattering rate 1 jr. Thus, one can 
expect that fluctuations of the Kondo temperature can exist 
even in the thermodynamic limit, when both A and E c van- 
ish. 

In order to explicitly connect the Kondo temperature to the 
LDOS fluctuations, we define as the Kondo temperature 
for a non-fluctuating LDOS, v, and use Eq. (01 into (Q~|i to 
find 2 ! 



T K = exp 



D ~ E * dE W^+M tanh 

w„ 2vE 



E 
(12) 



where Sp = p — v. From the fact that, on average, the devia- 
tions of the local density of states vanishes, (dp) = 0, and by 
expanding the right-hand side of Eq. (fT2l to second order in 
5p, we find 




D-E F rD-E P 

dE I dE' 

— Etp J — Err 



x ( tann | — — | tann 



2T 



(0) 

K , 



Sp(r,E F + E) 5p(v,E F + E'Y 



2vE 



2vE> 



(13) 



Note that Eq. ( fT3l provides an approximate expression for 
standard deviation of the Kondo temperature, 5Tk- When 



ST K < (T K ), we have ST K w T 



(0) 
K 



In" 



(Tk/T 



Using Eq. ©, the right-hand side of Eq. (\3[ can be broken 
down into three contributions, namely, 



T, 



K 



T 



(o) 



= 2 



K 



?,(r,A J Tf) + y,(r^ c ,rf) 
Q (t,E c ,T^]. (14) 



The function Sp arises from the spectral self-correlation term 
S(u>/ A) and is given by 



S,(r,A,T K ) = * (l + | 

D-E F 



1 + -Ren(0) 

P 



This term is of order A /Tr and therefore vanishes in the ther- 
modynamic limit as A — > 0. The decaying part of the spectral 
correlation function yields the second term on the right-hand 
side of Eq. ( TBi i. 



1 f rD—Ep dEdE i f E 



K 



xtanbYj-) [R 2 (E-E')-\] 



l + -R&TL(E-E') 



(16) 



In 2D, i?2 (w) — 1 decays as 1/uj for frequencies exceeding the 
Thouless energy. This causes Vp to be nonzero as E c — > 0, 
but this term turns out to be only of order 1 / g 2 . If correlations 
of wave functions are taken into account only up to first order 
in l/g, V/3 can be discarded. At finite E c , however, it needs 
to be taken into account since it yields terms of order E c /Tk 
due to spectral correlations at small frequencies, u> < E c . 
The term arising purely from the correlations of wave func- 
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tions is given by 

Q/3(t,E Ci T k ) 



1 

Aft 
x tanh 



d ~ Ef dEdE' , / E 
tanh 

/•;, \ - 1 h 

E' 



2T, 



K 



EE 1 "~""\2T R 
ReU(E-E'). (17) 



This term survives the thermodynamic limit. For instance, let 
us consider the 2D case. From Eq. (fTTT i. we find 



(18) 



For 1/t> Tk, one can use Eq. ( fT8l with E c = to write 
the double integral in Eq. ( TPTI i in terms of polylogarithm func- 
tions. To leading order and after setting Ep = D/2, we obtain 



Qf D) (r,0,T K ) 



Qngft 



In 



E F 

gTi 



A 



(19) 



which is valid when Ef/Tk 2> 3 3> 1. Combining Eqs. (fl4] i 
and ( fl"9l i, we arrive at 



hr 



X, 



A" 



T 



(0) 



2/) 



3n g [3 



In 



(20) 



which in entirely volume independent. 

Expressions for the (Tk) and STk can also be derived from 
Eq. ( fT3| ). However, in order to keep results accurate up to 
second order in Sp, one has to take into account the correlation 
between Tk and 5p when expanding the right-hand side of Eq. 
( fT2l . Doing so, we obtain for the average 



(T K ) = Tf {l 



1 



■ T^d T (a) 



Sp(r, 



A, T, 



(0) 



A" 



+Vp (r, £ c , T^) + Q,g (r, £ c , T^)] } . (21) 



Similarly, for the standard deviation, we arrive at 

2 r 



[5T K f = 2 (Tf ) (r,A,Tf ) +V ( 
+Q (r,£; c ,4 0) )" . 

Combining Eq. d20l ) with ( l22l we finally get 

r / \ -i 3/2 

1 / 1 \ 

$ t k\ 2 d 



t E T (0) 



(0) 



In 



(0) 



A 



(22) 



(23) 



This expression shows that even in the weak-disorder metal- 
lic regime (g 3> 1) the width of the distribution of Kondo 
temperatures remains finite as A and E c go to zero. It is in- 
teresting to note that this width is reduced by a factor 1 / ^/2 
when time-reversal symmetry is broken (j3 = 2). This effect 
can be understood if we recall that the weak-localization cor- 
rections, which account for exactly half of the wave function 
correlation correction, become suppressed when time-reversal 
symmetry is fully broken. Then, only wave function correla- 
tions due to the classical diffusion propagator remain. 



It is also important to remark that Eq. ( f23l is valid despite 
the fact that all electron states have a finite localization length 
£ in dimensions not exceeding two. In particular, we note that 
the width Eq. d23l exceeds previous estimates, such as 5Tk ~ 

J A(T { k\ which were based on the application of RMT on 
the scale of the local level spacing A^ and disregarded the 
dominant effect of wave function correlations that we consider 
here.^S 

From Eqs. d20b and OTT l. in the limit 1/t> Tk, the aver- 
age Kondo temperature is found to increase with the disorder 
strength as 



(Tk) 



2D 



T 



(0) 



A 



1 + 



1 



6Trg/3 



In 



1 



tT. 



(0) 
A , 



(24) 



Although the disorder-induced term in Eq. (124-b is only of 
order 1 /g, it is enhanced by the large factor coming from the 
third power of the logarithm of T^t. 

These calculations can be repeated using the diffusion prop- 
agator for quasi-one-dimensional (Q1D) wires. For open 
boundary conditions, the diffusion propagator in Eq. ( fTOb be- 
comes 



L/{k 2 F lA), 



oj<E c 



ReIl( Q1D >H = <^ V6/(k F A^¥[J), E c <u<1/t, 
{ 0, w > 1/r, 

(25) 

where L is the length and A is the cross sectional area of the 
quasi- ID wire. Therefore, in the long-wire limit, with E c <C 
Tk < 1 /V, we find 



In 2 



T, 



(0) 



A 



Q1D 



4tt\/3 1 1 



ft 



hp A 



tT 



(0) 



(26) 



A 



independently of the position of the Fermi energy. Hence, 
the average Kondo temperature increases with disorder by an 
inverse square root term, 



(Tj 



K/Q1D 



T 



(0) 
A 



ttV3 1 1 

~q~Wa 



tT 



(0) 



(27) 



A 



while the variance is nonvanishing by the square root of the 
same factor, 



ST K 



Q1D 



T, 



(0) 



A 



\ 



4ttV3 1 1 



kpA 



tT 



(0) 
A 



(28) 



The latter result is in agreement with the 1/g 1 ^ 4 dependence 
reported in Ref . 3 1 , where the asymptotic behavior for the 2D 
case, (J 3 /g) 1 / 2 , Eq. ( 1231 was also found. 



IV. NUMERICAL STUDIES OF A DISORDERED MODEL 

The analytical results presented in Sec. [HI] are applicable 
only in the weak-disorder regime and only predict the behav- 
ior of the two lowest moments of the distribution of Kondo 
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temperatures. In order to study more thoroughly the statisti- 
cal fluctuations of the Kondo temperature, as well as to check 
the accuracy of Eqs. (TfTT i. d22l i. d27l i. and (l28l l. we have car- 
ried out numerical studies of the fluctuations of Tk using a 
tight-binding model with nearest-neighbor hopping and ran- 
dom site potential to describe the conduction electrons. The 
model Hamiltonian reads 



N 



H = -t 



4 c j 



c\Ci 



(29) 



where each site potential Vi is drawn from a flat box distribu- 
tion of width W centered at zero. We assume each eigenstate 
of H to be spin degenerate, therefore sums over spins are im- 
plicit in Eq. (|29l . We consider only square lattices. Aspect 
ratios and boundary conditions are as follows: To simulate 
two-dimensional systems, we use a square geometry with pe- 
riodic boundary conditions along both directions; For quasi- 
one-dimensional systems, we adopt a rectangular geometry 
with hard-wall (periodic) boundary conditions along the short- 
est (longest) length. 

Employing standard numerical techniques, we have diag- 
onalized the Hamiltonian H for a large set of realizations 
of the disordered potential and different lattice geometries. 
The resulting eigenenergies E n and eigenvectors ip n (i), n — 
1, . . . , N, were used in conjunction with Eq. (Q~|i to determine 
Tk- In all simulations the Fermi level was placed at the lower 
quarter of the band in order to avoid the large peak in the den- 
sity of states at E = 0, reminiscent of the van Hove singular- 
ity found in the clean limit. No unfolding of energy levels was 
used. The simulations were carried out for 20 x 20, 10 x 100, 
and 8 x 200 lattice sizes, with 1 ,000 realizations for each value 
of W. In order to increase statistics, a total of 36 (20 x 20), 
196 (10 x 100), and 297 (8 x 200) different magnetic impurity 
sites were used. No magnetic flux was included, thus all re- 
sults refer to the time-reversal symmetric class. Partial results 
for the 20 x 20 case were previously reported in Ref. [l(| 

A connection between the tight-binding disorder strength 
W and the dimensionless parameter g (see Sec. [TTTb can be 
made by recalling the expression for the nonmagnetic scatter- 
ing rate in Born approximation, 



2ttv (V 2 ) 



(30) 



where v = 1/ AL 2 (we set the lattice constant a = 1, there- 
fore N = L 2 ). Noting that (V" 2 ) = W 2 /12 and A = D/N, 
we obtain 



1 _ 7T W 2 

t ~ 6~D~' 



(31) 
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Recalling that g = Ep t and setting Ep = D/A, we arrive at 
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FIG. 2: (Color on-line) Average (a), standard deviation (b), and their 
ratio (c) as a function of disorder strength for three values of J/D 
and three lattice geometries: two-dimensional (20 x 20) and quasi- 
one-dimensional wires (10 x 100 and 8 x 200). Dotted lines are only 
guides to the eyes. The dashed, dotted-dashed, and solid curves are 
the analytical prediction for J/D = 0.20, 0.30, and 0.50, respec- 
tively, based on Eqs. d21b and J22b and taking into account the finite 
values of E c and A. The analytical curves were plotted only within 
their region of validity. The same horizontal scale is used in (a), (b), 
and (c). 



Furthermore, when deriving Eq. QTl ) we assumed a flat den- 
sity of states (i.e., a parabolic dispersion relation), which is 
a rough approximation. Thus, taking into account such lim- 
itations, g has to be considered merely as a convenient way 
to parametrize the transition from a metallic (g ^> 1) to a 
strongly disordered, localized regime (g < 1). It is important 
to remark that g is not equal to the dimensionless conductance 
parameter^ = E c /A. 



D 

2tt \W 



(32) 



A. Average and standard deviation of Tk 



Equation (130b is valid only in the perturbative sense, when 
disorder is not too strong, and should break down near g = 1. 



The dependence on g of the average and standard deviation 
of the Kondo temperature, as well as their ratio, is shown in 
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Fig. 12 for several values of the exchange coupling J. The 
solid, dashed, and dotted-dashed curves are the analytical pre- 
dictions for a 20 x 20 (two-dimensional) lattice based on Eqs. 
(f2TT > and (l22l and taking into account the finite values of E c 
and A. The clean-limit value of the Kondo temperature, T^\ 
was obtained from Eq. ( IB II ) with the proper adjustment for 
the Fermi level position. 

Figure|2]contains several important features. First, the plots 
show that the statistics of Tk is quite independent of lat- 
tice size and geometry. Even wires with an aspect ratio of 
1:25 still present statistical fluctuations of Tk consistent with 
a two-dimensional geometry after the appropriate rescaling. 
Deviations are only apparent for the largest value of g con- 
sidered. This behavior can be understood if we recall that 
for diffusive dynamics the crossover from 2D to the Q1D 
regime occurs when the elastic mean free path becomes larger 
than the wire width but remains smaller than the wire length 
(L y < I < L x ). For a square lattice with weak site disorder 
(g ^s> 1), assuming a fiat density of states, it is straightforward 
to show that 




Applying this expression to the 10 x 100 and 8 x 200 lat- 
tices we find that the 2D-Q1D crossover happens around 
g — g c S3 7 and 6, respectively. Thus, at g < g c , the electronic 
eigenstates should follow the statistics of a two-dimensional 
disordered lattice, irrespective of the wire aspect ratio, and 
this is indeed what we observe. At g < 1, all states are local- 
ized and therefore the statistics of Tjc becomes independent 
of lattice size, geometry, or boundary conditions. 

In order to reach a regime where the analytical expres- 
sions for Q1D [Eqs. (l27T i and ( f28b l are applicable, we would 
need narrower wires or much weaker disorder. At the same 
time, we would need to maintain a large number of transver- 
sal modes and longitudinal diffusion. In practice, these con- 
straints can only be satisfied for much larger lattices than those 
we have investigated. 

It is also clear from Fig. |2]that the average Kondo temper- 
ature is only substantially affected by disorder when localiza- 
tion sets in (g — > 1). This behavior is reminiscent of the inde- 
pendence of the critical temperature on nonmagnetic disorder 
in weakly disordered conventional superconductors. 32 In the 
weak-disorder regime, (Tk) follows closely the bulk clean 
case value of the Kondo temperature [see Eq. ( IB H I. which 
is consistent with the prediction of Eq. (fJTJ. Deviations only 
become large in the strong disordered regime, where the ana- 
lytical expressions are not expected to be valid. The situation 
is less satisfactory for 5Tk, where the deviations from the an- 
alytical curves are seen for all disorder strengths. Note that 
the agreement with Eq. (l22l is better for large values of J/ D. 

Overall, the analytical expressions do capture the qualita- 
tive aspects correctly. There is an enhancement of the Kondo 
temperature and an increase of its variance as the disorder 
strength increases and g decreases. According to the analyti- 
cal calculations, this is mainly due to the appearance of weak 
correlations of wave functions at different energies. 



B. The distribution of T K 

The effect of disorder is much more drastic for the distri- 
bution of Kondo temperatures as a whole. In Ref. 10, it was 
found that a bimodal structure appears even for weak disorder 
and large exchange coupling. To further investigate this effect 
and its dependence on sample geometry, we present in Fig. [3] 
the distribution of Tk obtained from numerical calculations 
involving 10 x 100 and 8 x 200 lattices. These should be 
compared with Fig. 2 in Ref. [Tol where the 20 x 20 case was 
considered. 

Two distinct peaks are clearly visible when J/D is not too 
small. As J/ D increases, weight is transfered from the low- 
Tk to the high-TV peak. The distributions never come as 
close to a Gaussian as they do in the RMT case (see Ref. [13 
for the RMT distributions). The distinct features of the dis- 
tributions are a wide plateau in P(Tk) for weak disorder and 
a marked second peak at small Tk for strong disorder. Since 
Eq. ( TT9l > indicates that the enhancement of fluctuations in the 
Kondo temperature survive the thermodynamic limit, we be- 
lieve that the overall form observed in Fig. [3]should persist as 
the linear size of the lattice is increased. However, to verify 
this conjecture, one would need to perform a finite-size scal- 
ing analysis over a much wider range of lattice sizes. For an 
analytical proof, at least a few moments higher than the vari- 
ance would have to be calculated. Both are beyond the scope 
of this paper. The finite-size rescaling will be investigated in 
a future work using a more refined numerical technique^ 

It is important to remark that a finite-size scaling analysis 
was performed in Ref. 29 in the strong disorder limit only. It 
was found that a critical feature seen in Figs. Oc)-(f), namely, 
the divergence of the distribution as Tk — > at strong dis- 
order, does survive finite-size scaling. In Ref. 29, this diver- 
gence was found numerically to be a power law and it was 
argued that the exponent depended solely on the lattice space 
dimension. However, no analytical derivation was given to 
justify such a behavior. This issue is addressed in Sec. [V] 



V. DISTRIBUTION OF T K IN THE LOCALIZED REGIME 

For the purpose of investigating the power-law behavior re- 
ported in Ref. 29, let us consider the localized regime in 
more detail. When the localization lengths £ n of all states 
within the energy band are smaller than the system size, one 
can calculate the distribution function analytically by noting 
that only states which are located within a distance smaller 
than to the magnetic impurity state contribute significantly 
to the Kondo screening. The probability density of the nth 
eigenstate at the position of the magnetic impurity, ro = 0, 
is then given by (L/^ n ) d exp(— 2|r|/£ n ). We can then aver- 
age over the center-of-mass positions r„ of the states closest 
in energy to the Fermi level. Note that localized states can 
overlap. When the localization length is sufficiently large, 
several states can be found within a localization volume £ d 
and the corresponding energy levels will obey Wigner-Dyson- 
statistics. In particular, the two energy levels closest to the 
Fermi energy repel each other on the scale A^ = 
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FIG. 3: (Color on-line) The distribution of the Kondo temperatures for 10 x 100 and 8 x 200 lattices for a range of disorder strengths and 
three different values of J/D. Events corresponding to Tk = fall outside the vertical scale in (a) and (b). 



Under these assumptions, we can use Eq. (0) to derive a distri- 
bution of Kondo temperature in the localized regime, namely, 



P(T K ) 



Tk~ 



(F(e/T K )) ( 



(33) 



where is the sample volume, 

F(x) = 



tanh(z/2)cosfr(x/2) 



In 



d-l 



xT K 



D 



A ? tanh(a;/2) J 



(34) 

and (. . .) e denotes the average over the level spacing e be- 
tween the states within the localization volume £ d which are 
close to the Fermi energy Ep. Let us now look at both "weak" 
and "strong" localization regimes. 

For large localization lengths, is much smaller than the 
bandwidth and we can assume that energy levels in the vicin- 
ity of the magnetic impurity repel within the energy window 
A^. As a result, for Kondo temperatures exceeding the local 
level spacing, Tk > A^, the distribution becomes 



P(Tk) 



T K >At 



Tk 



In 



</-i 



Tk 
A? 



(35) 



which indeed contains a power law divergence at small Tk but 
with an exponent independent of the space dimension. As Tk 



becomes smaller than A^, this divergence is cutoff due to level 
repulsion. For a fixed level spacing e = A^, one gets an ex- 
ponential cutoff P{T K ) \T K< At~ (Af/Tl) exp(-A ? /T K ). 
However, the spacing e is known to follow the Wigner sur- 
mise, P(s = e/A^) = ap s" exp(— cp s 2 ), where, ap = n/2 
(32/tt) and cp = tt/4 (4/tt), for /3 = 1 (2). Performing the 
average over e, we arrive at 



K) \T K <A, 



Q ap Tj 



K 



(36) 



This distribution decays (rather than diverge) for Tk < A^. 

In the strongly localized regime, g -C 1, £ becomes of the 
order of the Fermi wave length Af- Then, only a single energy 
level has appreciable overlap with the magnetic impurity. The 
spacing between the energy level and the Fermi energy obeys 
a Poissonian distribution, P(e) = exp(— |e|/Af). Averaging 
over e yields 



P{T 



K) \T K >A e 



-Lin' 



d-l 



T 



K 



iF (37) 



and 



P{T K ) 



\T K <A E 



^/Tj 



(2 + As/T K ) 



In' 



d-i 



(38) 
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Thus, we also find a power law increase of P(Tk) with de- 
creasing Tk in the strong-disorder limit. But this divergence 
is also cutoff at Tk < Ag and P(Tk) converges to a finite 
value: P{T K ) \t k «a^ 1/A s . 

In summary, in contrast with the interpretation of the nu- 
merical scaling results in Ref. (I29I) . we find that the power 
of the divergence of P(Tk) does not depend on space dimen- 
sion, but is rather equal to 1. However, in both cases, the 
divergence is cut off at Tk < A^: For £ 3> A_f, there is a 
power-law decay, P(Tk) ~ T K , while for £ — * \p it ap- 
proaches a constant value, P(Tk) \tk<Ac~* VAg- These 
analytical results are found to be in qualitative agreement with 
our numerical results [see Fig. 0. 

It should be noted that while the events corresponding to 
Tk = are not always shown in Fig. [3] they were taken into 
account in the evaluation of (Tk) and STk in Figs. [2] The 
long tail toward the small-T^ region indicates that disorder 
enhances the probability of having unscreened, free magnetic 
moments at zero temperature. This is in agreement with the 
data present in Fig. 4 of Ref. where the probability of 
finding no solution to Eq. (Q~|i was plotted as a function of 
J/D. 



VI. DISORDERED GRAINS: RANDOM MATRIX THEORY 

In a small, weakly disordered, phase-coherent metallic 
sample, such as a metal grain or nanoparticle, 33 the single- 
particle energy levels repel each other and the wave function 
intensities are distributed randomly. This behavior becomes 
very relevant to several physical properties of the system at 
low temperatures, namely, when T < E c (see Fig. [TJ. This in 
turn will affect the Kondo temperature associate to magnetic 
impurities located inside the sample. 

For weak disorder, the dimensionless conductance param- 
eter is very large: Q — E c /A — ► oo. In this regime, the 
spacing between consecutive levels within a window of en- 
ergy ui < E c obeys the Wigner-Dyson distribution of RMT. 34 
Moreover, within a given symmetry class, the wave function 
intensities x n and the rescaled eigenenergies s n fluctuate in- 
dependently. The quantities x n , n — 1, . . . , N, are them- 
selves uncorrected and obey the so-called Porter-Thomas 
distributionsi22i2£ For the Gaussian unitary ensemble (GUE 
- broken time-reversal symmetry class) and for the Gaussian 
orthogonal ensemble (GOE - time-reversal symmetric class) 
these distributions are given by 

f p~ Xn GOF 
P ^ = { e-^/V2^ n , GUE. ^ 

In order to gain insight about the statistical properties of 
Tk over a wide range of values for J, we have solved Eq. 
(|2} numerically using eigenenergies s n obtained by diagonal- 
izing GOE and GUE random matrices of size N = 500. The 
resulting semi-circular spectrum was unfolded into a flat band 
following standard procedures and the Fermi energy was set to 
the middle of the band (Ep — 0). Instead of using the random 
matrix eigenfunctions to get the local wave function intensi- 




FIG. 4: (Color on-line) The dependence of average Kondo temper- 
ature on the exchange coupling for the GOE and GUE compared 
with the analytical solution in the clean case. Nonzero values of 
Tk are found below the clean-case threshold given by Eq. dB2l > 
(J- « 0.134 here). In this region, Tk is controlled by a few lev- 
els close to the Fermi energy, as a comparison with the results of the 
2LS model, Eqs. d45t and J46b . indicates. Inset: The dependence of 
the Kondo temperature average {Tk ) and standard deviation STk on 
the exchange coupling J for unfolded GOE spectra as compared with 
those obtained for a spectrum of equally spaced energy levels (picket 
fence). In both cases wave function intensities were drawn from the 
GOE Porter-Thomas distribution. It is only for (Tk) < A that the 
fluctuations of the energy levels result in an appreciable enhancement 
of (T K ) and ST k . 

ties x n , we generated these quantities directly from the Porter- 
Thomas distributions as given in Eq. (f39b- 35 A total of 500 
random matrix realizations were used for each ensemble type. 
For each realization (energy spectrum), we simulated differ- 
ent impurity locations by drawing 250 values of x n , thereby 
increasing substantially the statistics of Tjc without having to 
perform a large number of matrix diagonalizations. In Ref. 
WXk the distribution of Kondo temperatures obtained in this 
way were shown for several values of J/D. Here we provide 
a detailed analysis of these and other numerical and analytical 
results. 

In Fig. |U the GOE and GUE values of (T K ) are plotted 
versus J together with the curves obtained in the clean, ther- 
modynamic limit, namely, Eqs. ( 1BU and ( |B3t . The average 
Kondo temperature is found to coincide with the clean case 
regardless of the ensemble symmetry when (Tk) > A. It is 
only for (Tk) < A that it becomes dependent on the symme- 
try class. Note that spectral fluctuations allows for solutions 
of Eq. (0 below the clean-limit threshold value J_ (see Ap- 
pendixlB]). We will return to this point in Sec. I VI Bl 

A. Fluctuations in the large-Tx limit (Tk > A) 

The average and the standard deviation of the Kondo tem- 
peratures for GOE are shown in the inset of Fig. [4] They are 
compared with results obtained for an equally spaced spec- 
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trum (picket fence) where the eigenfunction intensities were 
also drawn from the GOE Porter-Thomas distribution. It is 
clear that the average and standard deviation of the Kondo 
temperature are insensitive to level fluctuations when Tk is 
larger than A. We can use this fact to obtain an approxi- 
mate expression for the distribution of Tk in the RMT: Using 
Eq. ( f39b and substituting the delta function by its Fourier inte- 
gral representation, we can perform the average over the wave 
function amplitudes in Eq. (0 exactly. For the orthogonal 
class, assuming an equally spaced spectrum, we obtain 



Pgoe(T k ) 
where, 



d \ ( DG r \ i r G v 

r—l x ' n—l.Tiy^r 



G n = (n- l/2)coth 



n - 1/2 
2k 



(40) 



(41) 



Equation (l40l can only be further simplified in some limiting 
cases, as outline below. 

For Tk ^S> A, one can set G r = r — 1/2 for r > K and 
G r = for r < k. The Kondo distribution then gains a 
stretched exponential form^ 



Pp (Tk) ~ 4- exp 
_ ™(o) 



(3k In ( — - ) — (3kq 



(42) 



where k = Tg> /A. Note that in the vicinity of T^> the 

distribution is close to a Gaussian: Near T^' and for «o > 1 
it can be approximated as 



Pp (Tk) ~ ~ exp [-/3(k - k ) 2 /2k ] 



(43) 



The departure from the Gaussian behavior occurs at the tails 
of the distribution: For k « /to (k » kq) the curve de- 
fined by Eq. d42b runs below (above) a Gaussian. In Ref. 
[Tol it was noted that P(Tk) becomes wider than Gaussian for 
small exchange couplings. Even for larger J, when the aver- 
age (Tk) exceeds A, clear deviations from the Gaussian be- 
havior were found, with the distribution showing asymmetric 
non-Gaussian tails. There was good quantitative agreement 
between Eq. ( l42l and the numerical data where kq > 1 and 
J < D. 

Using the Gaussian approximation of Eq. (l43l , we can 
establish a relation between the standard deviation and the 
average of the Kondo temperature in the limit Tk 3> A, 
namely/ ' 28 ! 36 



ST. 



K 1/3- 



{Tk) A 



(44) 



where (Tk) ~ T^ is assumed. Since (Tk) is ensemble in- 
dependent in the Tk S> A limit, one finds from Eq. (l44b that 
STk goe / STk gue = l/v2- These results are compared to 
the numerical data in Fig. [5] In general, when time-reversal 
symmetry is present ((3 = 1), level repulsion is weaker and the 
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FIG. 5: (Color on-line) The ratio ST K /(A{T K }) 1/2 as a function of 
the exchange coupling constant J for the GOE and GUE ensembles. 
Inset: The ratio between GOE and GUE values for (Tk) and STk- 
Analytical considerations predict these ratios to yield the universal 
numbers 1 and \/~2, respectively, over a wide range of exchange cou- 
plings J > J_. Deviations from this behavior are only seen for 
small exchange couplings, when J < J_ . These can be attributed to 
fluctuations of the two levels closest to the Fermi energy. 



tendency to localization stronger. Consequently, the probabil- 
ity of having a vanishing wave function at the magnetic im- 
purity position is enhanced, as well as the probability of large 
wave function splashes. This indicates that the RMT distribu- 
tion of Tk should be wider in this case than for the unitary 
class (broken time-reversal symmetry), in agreement with the 
results of Ref. The inset of Fig. [5] shows that, where 
applicable, the analytical prediction for the ratio between the 
standard deviations of the Kondo temperature for the orthogo- 
nal and unitary ensembles is in quite good agreement with the 
numerical results. For both RMT ensembles the distribution 
width scales with yA and therefore vanishes in the infinite 
volume limit. In Sec. |lll]we showed that this is no longer true 
when spectral correlations beyond the Thouless energy scale 
exist. 

There are realizations where no solution with Tk > can 
be found for Eq. (0. They can be interpreted as the absence 
of screening of the magnetic impurity dipole moment. For a 
clean system, one can estimate the minimum exchange cou- 
pling J_ where Eq. (O ceases to provide a finite value for 
Tk (see Appendix [B). Non-magnetic disorder strongly af- 
fects J_ . For small exchange couplings, free moment events 
amount to a large portion of all events. Their probability is 
plotted separately in Fig. [6] It is remarkable that even for 
J < J_, where (Tk) <C A, the distribution shows a clear 
maximum at a finite value of Tk- The distribution (see Ref. 
1 1 QD decays towards Tk — even though there is a large 
amount of free moments in this case (see Fig. [6]). This peculiar 
behavior can be qualitatively understood by recalling the level 
repulsion of RMT ensembles and can be quantitatively treated 
within a two-level model, as outlined in Sec. IVIBI below. 
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FIG. 6: (Color on-line) Probability of finding free magnetic mo- 
ments. The data points correspond to the fraction of events in the 
RMT simulations where either Tk = or no solution to Eq. 10 was 
found, and is plotted as a function of the exchange coupling J. The 
dotted line indicates the point where J = J_ (clean case threshold). 
There is a good agreement between simulations and the analytical 
predictions for J< J- based on the two-level system model. The 
inset shows the tail of the probability function at J > J_ . 



system for a given bare exchange coupling. This function is 
equal to the probability that the Kondo temperature is either 
zero or that there is no solution to Eq. (0. In these cases, the 
correction to the exchange coupling remains small for all tem- 
peratures and the Kondo effect is quenched, in the sense that 
the magnetic moment remains unscreened down to zero tem- 
perature. As discussed in Appendix|B] in the clean limit the 
free moment probability is a sharp step function of J, namely, 
-Pfrcc(</ <•/-)= 1, and Pf ree (J > J_) = 0, with J_ given 
by Eq. dB2b . When there is randomness, J_ is a function of 
the random wave function intensities x n and level spacings s n 
of the closest levels to the Fermi energy. Upon averaging over 
their distributions, we expect the step function to become a 
smoothly decaying function. Thus, due to the finite probabil- 
ity that level spacings at the Fermi energy exceed A or wave 
function intensities are small, there is a finite probability of 
finding free moments even at J > J_, where the magnetic 
moment would be screened without randomness. We have de- 
termined this probability distribution numerically (see Fig. [6]). 
For J < J_, we can compare it with the result obtained within 
a 2LS approximation. A first attempt can be made by setting 
the wave function amplitudes to be constant, xi = 1 and fix- 
ing all but the two closest levels to the Fermi energy, which 
are allowed to fluctuate according to the Wigner surmise. 34 
For the GOE we obtain 



B. Two-Level Approximation (Tk < A) 

For Tjc < A, only the two states closest in energy to the 
Fermi energy are appreciably coupled to the magnetic impu- 
rity. For small exchange couplings J < J_, this affects the 
whole distribution of Kondo temperatures, while for larger ex- 
change couplings it applies only to the small-Tjf tail. In this 
limit, we can simplify the analytical calculations by separat- 
ing the contributions to the sum in Eq. (O into two groups. 
In the first we include only the two nearest levels to the Fermi 
energy and fully take into account the fluctuation of their rela- 
tive spacing and their wave function amplitudes. The remain- 
ing levels form the second group and are treated assuming 
an equally spaced spectrum and x n = 1. When the wave 
functions at the Fermi level are distributed according to the 
GOE Porter-Thomas distribution but the energy levels are kept 
fixed, the average Kondo temperature is found to be 



-Ei(-l/2) / D D 

(Tk)goe = Ye Aexp - 7 + — 



(45) 



where Ei(x) denotes the exponential integral function,— with 
Ei(-l/2) w -0.56. For the GUE Porter-Thomas distribu- 
tion, we obtain instead 



(Tk. 



GUE — 



-2Ei(-l/2) / 2D 



J 



2D 

1Z 



(46) 



showing that the average Kondo temperature in the unitary 
ensemble becomes exponentially smaller than in the orthog- 
onal ensemble, which is consistent with the numerical result 
presented in Fig. [5] 

It is of particular interest to determine the probability 
Pfree(^) of having unscreened, free magnetic moments in the 



Pfree(J) = exp 



(D/J-D/J- +2) 2 



(47) 



which is valid only for J < J + = J_/(l — J-/D), since for 
larger values of J the fluctuations of the other levels can no 
longer be neglected. 

When the wave function amplitudes are constant, the prob- 
ability to have free magnetic moments is just proportional to 
the probability that the level spacing at the Fermi energy is of 
the order of (J). That is exponentially small for J > J_. 
However, in Fig. [6] we see that the decay with J is slower than 
that, indicating that the fluctuations of wave functions are cru- 
cial. Indeed, under realistic conditions, there is no appealing 
reason to take xi = 1. When wave functions are allowed 
to fluctuate according to the Porter- Thomas-distribution, we 
rather obtain (again for the GOE) 



P' 



GOE 



free 



(J) 



dx exp I —x — 



uj exp 



Erfc 



(48) 



with 



uj 



D D 

H 

2J 23- 



(49) 



Note that P£° E (0) = 1 and P£° E ( J J) « 2/tt. For J > J_, 
the fluctuations of the other energy levels become important 
and the 2LS approximation can no longer be used. Similarly, 
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for the GUE we obtain 



~ Jo 



UjX) 



x exp 



-UJX 

7T 



(50) 



These results are compared to the numerical calculations in 
Fig. [6] and good agreement is found for J < J_. Thus, in 
RMT, we can conclude that fluctuations of wave functions do 
enhance the probability of finding free moments for J > J_ . 

Finally, there is remarkable feature of the distribution of the 
Kondo temperature, as shown in Ref. Oil especially notice- 
able for small exchange couplings (J/D < 0.1): Although the 
number of free moments is large in this regime, the distribu- 
tion still has a distinct maximum and turns to zero as Tk — > 0. 
We find that this is due to the energy level repulsion of RMT. 
Evaluating the distribution function in the 2LS approxima- 
tion (i.e., for a fixed energy level spacing but fluctuating wave 
functions), we arrive at P{T K ) ~ (A/T K ) e - A / 2TK . Then, 
using the Wigner surmise to average this expression over level 
spacings, we find that the distribution decays as a power law, 



P(T K < A) 



T, 



K 



(3 + 2A^ 



(51) 



with ax = 7r/2 and a 2 = 32 /tt, in fair qualitative agreement 
with the numerical data presented in Ref. 



VII. DEPENDENCE OF P{Tk) ON THE 
CONCENTRATION OF MAGNETIC IMPURITIES 

In the absence of an external magnetic field, the symme- 
try class of the underlying single-particle basis to the Kondo 
problem is controlled by the concentration of magnetic im- 
purities and the relation between certain energy scales. The 
general idea is that, as the number of magnetic impurities 
increases, statistical fluctuations of the electron states in the 
sample cross over from the orthogonal class (time-reversal 
symmetric) to the unitary class (broken time-reversal sym- 
metry). This is due to the fact that the spin dynamics of the 
magnetic impurities can be slow compared to the time scale 
of the conduction electrons, thereby breaking effectively the 
time-reversal invariance on their time scale^iS However, as 
we will argue below, the dynamical regime in the sample, and 
consequently the strength of the Kondo screening, also plays 
a crucial role in determining on which side of this crossover 
the system finds itself^ - 

We begin by recalling that for a mesoscopic sample in the 
weak localization (WL) regime, i.e., at energies E > E c 
(short time scales), the dimensionless parameter controlling 
the orthogonal-unitary crossover due to magnetic impurities 
can be written as^ 



X 



WL 



1 



(52) 



E c t s : 

with the crossover centered at X s = 1. As the energy is 
lowered below E c (long time scales), one enters the zero- 
dimensional, RMT regime, where the average level spacing 



A takes over E c as the relevant energy scale. In that regime, 
the crossover parameter is given by 



X 



RMT 



1 

a7 



= gx 



WL 



(53) 



instead. Thus, for a sample in the RMT regime, the crossover 
occurs for spin scattering rates smaller by a factor 1/Q with 
respect to the weak-localization regime, where Q is the dimen- 
sionless conductance. For the RMT regime, this crossover 
has been recently studied for a Kondo quantum dot in an 
Aharonov-Bohm ring^i 

At any finite temperature, the spin scattering rate l/r s is 
renormalized by Kondo correlations: It is small at both T ^S> 
Tk and T <C Tk, having a maximum at around T = Tk, 
as has been observed in weak-localization experiments! 4 ^ 3 . 
This behavior is analogous to that observed in the reentrance 
of gapless superconductivity. 44 For a clean sample (see Refs. 
30.31, 



Ts{T) 



7rn m 5(5+1) 
(TV 



In" 2 {T/T K ), T>T k , 

TkT k , ( (54 ) 



10 



-cfl, 



T^T K , 



where n m is the density of magnetic impurities, S denotes 
their spin, a is a numerical factor smaller than unit, and 
w « 0.41 is the Wilson number. In Refs. andl31l it was 
found numerically that a as 0.2. In Eq. d54l . cfl is a factor 
that arises from the relation between the inelastic scattering 
rate and the temperature-dependent dephasing rate^ In 2D, 
cfl ~ 0.946. 

A few elucidating remarks about Eq. d54l ) are in order. 
First, we define the spin scattering rate as in Ref. 38, which 
is smaller than the dephasing rate defined in Ref. [3l| by a 
factor 1/2. Secondly, note the difference between our defi- 
nition of Tk and that used in Ref. 3 1 : Ours corresponds to 
their perturbative Tk- Thirdly, the numerical prefactor of the 
(T/Tk) 2 term in Eq. ( l54b was obtained within Fermi liq- 
uid theory, 1 where the exact result differs from that obtained 
in the standard Sommerfeld low-temperature expansion by a 
factor 3. In the following, we will approximate the numer- 
ical prefactor multiplying the (T /Tk) 2 term by unit, since 
w 2 7t 4 c F l/16 = 0.968. 

The low-temperature limit of the spin scattering rate shows 
the expected Fermi-liquid scaling based on Nozieres' the- 
ory for the Kondo problem, which is valid at T <C Tk*& 
At temperatures exceeding Tk, Eq. ( 1541 is consistent with 
the perturbative poor man's scaling i 2 ' 21 ' 44 ' 46 ' 47 However, we 
note that recent experiments have shown the scattering rate to 
be approximately linear with temperature for a wide interval 
below Ta'. 48 This behavior has been explained theoretically 
by Zarand and coworkers using a numerical renormalization 
group calculation of the frequency-dependent inelastic scat- 
tering rate, 30 where they also obtained a w 0.2. 

In principle, one could expect that the maximum value 1 /r s 
can reach should be given by the unitary limit of the scatter- 
ing cross section. However, in reality, the maximum value is 
found to be smaller than that by the factor a shown in Eq. 
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(O. In 2D, 30 



(55) 



Thus, the maximal value the crossover parameter can reach 
for a two-dimensional sample in the weak localization regime 
is 



X 



WL 



N„ 



= a — 

max 



(56) 



where N m = n m L 2 is the number of magnetic impurities for 
a sample with linear size L. When there are only a few mag- 
netic impurities, N m < Q, the crossover parameter is small 
and the sample is in the orthogonal regime. Increasing the 
concentration of magnetic impurities increases the parameter 
X s and eventually leads to the unitary regime. As found in 
Sec. I VII this is accompanied by a decrease in the width of the 
distribution of Kondo temperatures. 

It has recently been pointed out that the correlations be- 
tween wave functions at different energies are enhanced in 
the GOE-GUE crossover regime. 16 Thus, since the width of 
the distribution of the Kondo temperature is enhanced by wave 
function correlations, one can expect a widening of the Kondo 
distribution in the GOE-GUE crossover regime due to the 
presence of a weak magnetic field or a small amount of spin 
scattering. 49 The wave function correlation function is then 
given by 



2A 2 A 2 



=E n -E„ 



4A 4 A 2 + tt 2 lu 2 ' 



n 2 (\Mr) \ 2 \^m(r) | 2 )| 

(57) 

Here, A is a dimensionless parameter related to the magnitude 
of the time-reversal breaking perturbation. In the presence of 
an external magnetic field B, X 2 = l/r^A, where 1/tb is 
the magnetic phase-shift rate. For a diffusive quantum dot, 
this rate is given by 1/tb = je 2 D e B 2 L 2 /h 2 , where 7 is 
a geometrical factor of order unity. In the presence of mag- 
netic impurities, the GOE-GUE crossover is governed by the 
spin scattering rate and A 2 = l/r s A. Remarkably, the mix- 
ing of GOE eigenfunctions brought by time-reversal symme- 
try breaking induces wave function correlations over a macro- 
scopic energy scale equal to either 1/tb or l/r s . However, 
the amplitude of these correlations is proportional to t s A and 
therefore they vanish in the thermodynamic limit. 16 In order to 
see that, we use Eq. ( PT21 to explicitly evaluate the width of the 
distribution of the Kondo temperature related to the crossover 
wave function correlations. We obtain 



(5T 2 k) 



1/(1-1 



dt , 1 + 1 

— In . 

t 1-t 



(58) 



Thus, for > l/r s , these GOE-GUE correlations only 

contribute to a width 8Tr ~ a/t^A/V that vanishes in 
the thermodynamic limit when A — > 0, but is larger than 
the value obtained in the pure RMT ensembles by a factor 



As the magnetic impurity concentration increases further, 
the superexchange interaction between magnetic impurities 
begins to compete with the Kondo screening. 50 The superex- 
change interaction coupling Jrkky fluctuates and its aver- 
age and spreading depend on the amount of disorder, impu- 
rity concentration, and details of the Fermi surface. Jrkky 
is zero on average, 51 but fluctuates according to a wide log- 
normal distribution. 52 Its typical value is of the same or- 
der as that of a clean sampled namely, \/ (Jrkky) = 
(n m / v)(J / D) 2 cos(2kpR). Even when the typical superex- 
change coupling constant is smaller than Tk, there is a small 
chance that clusters of localized spins form. When two lo- 
calized spins couple ferromagnetically to form a triplet state, 
they contribute to the dephasing rate of itinerant electrons. 
Such a contribution to the dephasing rate scales with n^M^ 
When Jrkky exceeds the Kondo temperature Tk of a sin- 
gle magnetic impurity, the spins of the magnetic impurities 
are quenched. They form a classical random spin array whose 
spin scattering rate is smaller than that for free quantum spins 
by the factor 1/3, and scale linearly with n m . 55 At interme- 
diate concentrations of magnetic impurities Griffiths-McCoy 
singularities can appear and induce an interesting though more 
complex behavior (see Ref. 56 for a review). 



VIII. DEPHASING DUE TO FREE MAGNETIC MOMENTS 
IN DISORDERED METALS 

We now discuss the relevance of the distribution of the 
Kondo temperature Tk to the quantum corrections of the con- 
ductance of mesoscopic wires. 

As discussed in Sec. IVHI [see Eq. d54l >, in particular], the 
temperature dependence of the spin scattering rate has a max- 
imum at Tk, resulting in a plateau of the dephasing time^Sd 3 . 
However, we have found in Sec. [IV] that disorder can result 
in a finite probability of having magnetic impurities with van- 
ishingly small Kondo temperatures. Thus, the question arises 
whether the latter effect can yield a finite contribution to the 
dephasing rate at temperatures below the average Kondo tem- 
perature. It has recently been argued that the dephasing rate 
due to a magnetic impurity can be related to the inelastic cross 
section of a magnetic impurity at the Fermi energy^^^ 



1 



n m VF&inc\, 



(59) 



where vf is the Fermi velocity, and the inelastic cross section 
can be obtained from the difference between the total and the 
elastic cross section, 



Cincl(-E^ 



Ctotal — Oel- 



(60) 



Using the optical theorem, both the total and the elastic cross 
sections can be related to matrix elements of the T-matrix of 
the magnetic impurity. 30 Based on the fact that in the deriva- 
tion of the weak-localization correction one averages over the 
disorder potential and thereby recovers translational invari- 
ance, the authors of Refs. |3(J and[3l| calculated the spin scat- 
tering rate of Eq. d59l for a plane waves basis. However, one 



should calculate the spin scattering rate as averaged over all 
momentum directions in order to correctly determine how the 
dephasing rate caused by magnetic impurities is affected by 
the random environment. Following this approach, we find 
that the spin scattering rate of a state of energy E should be 
given by 57 

l/r s (E, T) = ^J2 S ( E p ~ E ) [ Im 



(61) 



p' 



Using the relation between the T-matrix element and the 
propagator of a localized d-level in the Anderson model, 30 
G d (E n ,T), we obtain 



(p\T\p') = -Sli; p ,(0)r p (0)t 2 G d . 



(62) 



Inserting this relation into Eq. (|6TT >, we find that the spin scat- 
tering rate is given by 



l/r s {E,T) 



where 1 / Ts is the inelastic spin scattering rate in a clean sys- 
tem. 1/ts°\e,T) has a universal dependence on T/Tk, as 
has been shown in Refs. l3(il andl3ll Thus, Eq. d63l leads us 
to conclude that the spin scattering rate not only depends on 
the ratio T/Tk but also explicitly depends on the local den- 
sity of states. The same conclusion has been recently reached 
independently in Ref. (58). But both Tk and p(E, 0) are ran- 
domly distributed, taking different values for each magnetic 
impurity, and so will 1/t s (T). Note that the distribution of 
Kondo temperatures alone is not sufficient for determining the 
average spin scattering rate. We have seen in Sec. [Till that in 
a two-dimensional metal, the wave functions correlations at 
different energies are of order 1/g. The probability of having 
simultaneously a small Kondo temperature and a small local 
density of states at the Fermi energy leading to an anomalous 
value for the spin scattering rate is thus expected to be of order 
1/g as well. 

We can derive the temperature dependence of the spin scat- 
tering rate in a sample with a finite number of magnetic im- 
purities by averaging the dephasing rate of a single magnetic 
impurity over an ensemble of magnetic impurities. To facil- 
itate the calculations, we approximate the universal function 
for the dephasing rate by the following function, 



p(E,0). 



p(E,0) 



(4 lmG d + wtt G 2 d ) 



i/rj o) (£,n 



(63) 
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-ir 2 S(S+l) 
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, (64) 



with a — 0.2, as found numerically in Ref. |3CJ Note that tem- 
perature scales with Tk only. Equation ( f64b coincides with 
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FIG. 7: The spin relaxation rate as a function of temperature for 
J/D — 0.24. The broken lines correspond to ensemble averages of 
Eq. (63) over fluctuating Kondo temperatures and the LDOS at the 
Fermi energy (data from the 20 x 20 lattice simulations was used). 
The solid line is obtained from Eq. (64) setting Tk = {Tk ) (clean 
case). 



Eq. d54l > in both low- and high-temperature limits. In Fig. [7] 
we show the result of using Eq. d64l ) in Eq. d63l and ensem- 
ble averaging the spin relaxation rate over Tk and p(Ep, 0) 
obtained numerically for the 20 x 20 lattice model (see Sec. 

One finds a clear departure from the clean case prediction 
at low temperatures^ which amounts to an enhancement of 
the spin scattering due to the presence of nonmagnetic dis- 
order. Moreover, the maximum value of the spin scatter- 
ing rate is found to be reduced in proportion to the disorder 
strength and its position shifted to larger temperatures. A 
very similar behavior has been recently observed in weak- 
localization measurements performed by the Grenoble group 
on quasi-one-dimensional Ag wires doped with 2 to 20 ppm 
Fe impurities^ as well as by the Michigan group, 60 where 
similar samples where implanted with 2 to 10 ppm of Fe im- 
purities. One should note that these samples have very high 
diffusion constants, of order D e = 300 — 400 cm 2 /s, so that 
the parameter g in these samples is large, g > 100. For the 
Fermi energy in Ag, Ep — 5.5 eV, and the Kondo temper- 
ature of Tk = 4 K, the ratio between band width and ex- 
change coupling is found to be of the order of J/D = 0.1. 
Thus, according to Eq. ([23), ST K /{T K ) < 1 is fulfilled and 
the width of the distribution of Kondo temperatures should 
be negligible in 2D samples. However, one has to be careful 
when considering the relevance of our results to the weak- 
localization experiments of Refs. |H|and|6ll F° r instance, we 
note that the low-temperature dephasing rate is dominated by 
the non-Gaussian Iow-Tr- tail of the distribution of the Kondo- 
temperature, which we have found to be present even for weak 
disorder, g 3> 10. Moreover, the samples used in Refs. Inland 
60 were quasi-one-dimensional, in which case 5Tk is known 
to be further enhanced (see Sec. ITTTb - 

In order to connect theory with experiments, some more 
work needs to be done. On the theoretical side, a system- 
atic study of how the distribution of the Kondo temperature 
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scales with the width and length of the sample is necessary. 
This is the subject of ongoing work and will be published 
elsewhere. 11 On the experimental side, the relevance of the 
distribution of the Kondo temperature to the low-temperature 
anomaly of the dephasing rate could be established by utiliz- 
ing samples with lower diffusion constants, such as the low- 
mobility samples of AuPd, examined in Ref. 61 (see also Ref. 
l62t) . These have diffusion constants as low as D e = 13.4 
cm 2 /s, corresponding to g w 17. The presence of free mag- 
netic moments would be facilitated in semiconductors such as 
Si and GaAs where the effective mass of electrons is smaller 
than in metals by a factor 30. Moreover, 2D electron gases 
with akp ~ 1 can be produced with those materials. In 
low-mobility GaAs wires, Anderson localization has been ob- 
served even in wires with g 30. 63 Thus, all mesoscopic en- 
ergy scales can become relevant when the system reaches tem- 
peratures of the order of 1 K. However, so far no magnetic im- 
purities with detectable Kondo temperatures have been found 
in semiconductors, since the low density of states suppresses 
the Kondo temperature exponentially. 



IX. CONCLUSIONS 



low-dimensional disordered systems since the variance of Tk 
depends markedly on the disorder strength. This conclusion, 
however, is at odds with the proposal of Ref. |65[ where it was 
argued that in the low-temperarature, strong coupling limit, 
the Kondo effect is not affected by disorder at all. We note 
that this result was obtained using a path integral representa- 
tion and a semiclassical expansion (the large-A^ limit is taken, 
where N — 2s + 1 is the number of spin components). The 
disagreement between Ref. |65| and our results may not be 
real, however, since their calculation is limited to T <C Tk 
and cannot rule out the existence of magnetic impurities with 
exceedingly small Tk ■ 

It remains to be studied how the distribution of Tk scales 
with system size and how it changes (if at all) in the quasi-one- 
dimensional limit. The numerical techniques employed in this 
work did not permit the investigation of sufficiently large lat- 
tices. It also remains to be checked whether the large fluctua- 
tions of Tk seen at strong disorder survive a non-perturbative 
approach such as the numerical renormalization group. Fi- 
nally, we note that we have not attempted to explore the con- 
sequences of having a free magnetic moments to the nature of 
the electron liquid. This issue has been recently addressed by 
other authors in the context of heavy-fermion materials.^S 



The screening of magnetic moments in metals, the Kondo 
effect, is found to be quenched with a finite probability in the 
presence of nonmagnetic disorder. For weak disorder, g 1, 
the effect is shown analytically and numerically to be due 
to wave function correlations. Albeit weak in the metallic 
regime, these correlations have a sizeable effect on the Kondo 
temperature. That is because the Kondo temperature involves 
a summation over all eigenstates with a finite wave function 
amplitude at the location of the magnetic impurity. As dis- 
order increases and the sample goes from the diffusive to the 
localized regime, an increase in the correlations of the LDOS 
at different energies induces even stronger fluctuations in the 
Kondo temperature. Thus, the distribution of the Kondo tem- 
perature P(Tk) retains a finite width in the limit of vanish- 
ing level spacing A. When time-reversal symmetry is broken 
either by applying a magnetic field or by increasing the con- 
centration of magnetic impurities, P(Tk) becomes narrower. 
The probability that a magnetic moment remains free down 
to the lowest temperatures is found to increase with disorder 
strength. This result is consistent with direct measurements of 
Kondo effect in Cu(Fe) thin films and Au(Fe) wires.— 

We have also shown that magnetic impurities with a small 
Kondo temperature are shown to modify the temperature de- 
pendence of the dephasing rate at low temperatures, T <C 
(Tk), as measured in weak-localization transport experi- 
ments. 

Equations ( f20b and d26b allow us to conclude that the dis- 
order can substantially affect the Kondo temperature. This 
could be contrasted with the effect of disorder on another en- 
ergy scale derived from a many-body state, namely, the criti- 
cal temperature of a superconductor, T c . Anderson's theorem 
states that T c is not affected by nonmagnetic disorder.— Our 
results demonstrate that there is no analog to this theorem for 
the Kondo temperature of magnetic impurities embedded in 
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APPENDIX A: DERIVATION OF EQ. Q 

We begin with the Anderson model, namely, a localized d- 
level coupled to a conduction band^L 

H = ^ E nn n <r + £d ^ " dCT + U ^d+fld- 
n,a a 

+ ^ {tndCnaCda + tdnC^ a C n< y) , (Al) 

71, (X 

where E n are the eigenenergies of the conduction band elec- 
trons, Ed is the energy of a singly occupied d-level, and Sd + U 
is the energy of a doubly occupied d-level. td n = t* d are 
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the hybridization matrix elements between the d-level and the 
conduction band state Projecting out of the Hilbert space 
all state where the <i-levels are doubly occupied one obtains 
the Kondo Hamiltonian 



Hj = Js- S, 



(A2) 



where the matrix elements of the exchange interaction in the 
basis of the eigenstates of the conduction electrons are given 
by 



Jkl = 



8 tkd tm 
U ' 



(A3) 



These matrix elements are positive and thus the interaction is 
antiferromagnetic. The hopping matrix element connecting 
the localized rf-state 0d(r) to the conduction band state ^> n (r) 
is given by 



d d rr d (r)V(v)Mr), 



(A4) 



where V is the potential energy of the tunneling barrier. For 
an impurity state strongly localized at r = we obtain 



J u « J^(0)V«(0), 



(A5) 



where J = 8(1 / 'm* o,q) 2 U ', m* is the band mass, and ao is the 
radius of the localized state at the magnetic atom. We define 
Tk by the divergence of second-order perturbation theory. To 
second order in J, there are two processes to be considered: 
(i) The scattering due to the exchange coupling J of an elec- 
tron from state |n) to a state |Z) close to the Fermi energy via 
an intermediate state |m). This process is proportional to the 
probability that state \m) is not occupied, 1 — f(E m ), where 
f(E) is the Fermi distribution, (ii) The reverse process, in 
which a hole is scattered from the state \l) to the state |n) via 
the occupied of a state |m) with probability proportional to 
the occupation factor f(E m ). Thus, we find that the exchange 
coupling is renormalized to 
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2N E„ 



r»|^ m (o)| 2 



E f 



tanh 



E m — Ef 

2T 



(A6) 

For positive exchange coupling, J > 0, perturbation theory 
diverges as the temperature is lowered. Defining the Kondo 
temperature as the temperature where the second-order cor- 
rection to the exchange coupling becomes equal to the bare 
coupling, we arrive at Eq. (|T). 

An equivalent expression can be derived from the renormal- 
ization group equation, which in the two-loop approximation 
is given byS 



dJ 
dt 



-r- = J 



2 p{E F + A)+p{E F ~ A) 
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2D 2 v 
1 J 3 p(E F 
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A dE p(E F + E)/v 
„ X (1 + E/A) 2 

A dE p(E F - E)/v 
„ X (1 + E/A) 2 ' 



(A7) 



where t = lnD/ (2A). Defining Tk as the value of A at which 
J flows to the band width D and solving Eq. dA7| ) for Tk, we 
obtain for a clean system 



r (0) 

[ _fC2-loop 



e J 



(A8) 



However, if we keep only the first term on the right-hand side 
of Eq. ( IA7t - the one-loop term, we obtain a self-consistency 
equation for the Kondo temperature that reads 
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T; 



K 



T, 



(0) 
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dt 



2 ■'•ZT K /D<\t\<l t 



5p{E F + Dt/2) 



(A9) 



For Tk ^> A, Eq. JA91 > coincides with the self-consistency 
equation obtained using perturbation theory, namely, Eq. (fl]i. 



APPENDIX B: CLEAN CASE 

We briefly describe the behavior of the Kondo tempera- 
ture in the clean, bulk limit. For that purpose, we assume a 
spectrum of N equally spaced levels, band width D = NA, 
and spatially uniform wave function intensities (plane waves). 
When the Fermi energy is in the middle of the band and is 
even, all levels are either doubly occupied or empty at T = 
and we have s n =n- N/2 - 1/2, n = 1, . . . , N. For N > 1 
and Tk ^> A, we then find for Eq. (fTJ the well-known solu- 
tion 



T^ 0) w 0.57L>exp(-L>/J), 



(Bl) 



which agrees to lowest order in J/D with results from 
more accurate methods such as the numerical renormalization 
group. In fact, for the latter, the next leading order correction 
in the exponent has been found to be — . 5 ln(D/ J) 4 indicat- 
ing that our treatment is valid for J < D up to preexponen- 
tial corrections. For small J, however, Tk approaches A and 
turns abruptly to zero at 



J_ 



D 



ln(27Y) + C 



(B2) 



in a nonanalytical fashion, namely, 



Tf (J- J_) 



A 



2 ln[(D/J_ - D/J)/A] 



(B3) 



(here C w 0.58 is the Euler number). For J < J_, Eq. (fl} 
has no solution in the clean limit. That means that the lowest 
order corrections to the bare exchange coupling remain weak 
down to zero temperature and marginal terms have to be in- 
cluded in order to determine Tk- 4 Nevertheless, Tk <C A in 
this regime and we can conclude that there will be unscreened 
magnetic impurities (free moments) down to the lowest acces- 
sible temperatures in bulk clean metals when J < J_ . 
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